Nonlinear resonance reflection from and transmission through a 
dense glassy system built up of oriented linear Frenkel chains: 

Two-level model 

O 
O 

<N 

5_! 1 E. Conejero Jarque 

q.: 

Departamento de Ftsica Aplicada, Universidad de Salamanca, E- 37008 Salamanca, Spain 

m 

V. A. Malyshev 

National Research Center "Vavilov State Optical Institute", Birzhevaya Liniya 12, 199034 
' Saint-Petersburg, Russia 

O 
O 

• i-H 

C/> 

>Y 

, Ch'. Abstract 



Oh- 



(February 2, 2008) 



' A theoretical study of the resonance optical response of assemblies of ori- 

<^ ■ ented short (as compared to an optical wavelength) linear Frenkel chains is 



carried out. Despite the fact that the energy spectrum of a single chain is 
composed of the bands of Frenkel exciton states, a two-level model is used 
to describe the optical response of a single linear chain. We account for 



only the (on-resonance) optical transition between the ground state and the 



^ ■ state of the one-exciton band bottom as having the dominating oscillator 

5-H ' 



strength as compared to the other states of the one-exciton manifold. The 
(off-resonance) process of creation of two excitons per chain is neglected be- 
cause it requires a higher excitation frequency due to the quasi-fermionic 
nature of one-dimensional Frenkel excitons. A distribution of linear chains 
over length resulting in fluctuations of all exciton optical parameters, such as 
the transition frequency and dipole moment as well as the radiative rate, are 
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taken explicitly into account. 

We show that both transmittivity and reflectivity of the film may behave 

in a bistable fashion, originated from saturation of the nonlinear refraction 

index, and analyze how the effects found depend on the film thickness and on 

the inhomogeneous width of the exciton optical transition. Estimates of the 

driving parameters show that films of oriented J-aggregates of polymethine 

dyes at low temperatures seem to be suitable species for the experimental 

verification of the behavior found. 
PACS number (s): 42.65.Pc, 78.66.-w 
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I. INTRODUCTION 

Since pioneering works of Jelley§ and ScheibJ, who discovered the phenomenon of linear 
aggregation of polymethine dye molecules in solution with rising their volume concentration 
and, as a result, drastic changes of their optical properties (the appearance of a red-shifted 
narrow band, called now J-band), these species have been the subject of unremitting interest 
of reseachers. In the beginning of the nineties, Wiersma and co-workersiHl demonstrated the 
possibility of further narrowing of the J-band as well as shortening of its emission lifetime 
by an order of magnitude via cooling to the temperature of liquid helium. At present, it 
is widely accepted that the unusual properties of such systems originate from the fact that 
their optically active states are Frenkel exciton states (see for a comprehensive review Refs. 
|and 0). 

In recent time, a considerable attention has been drawn to the problem of strong cou- 
pling of organic polymers and molecular aggregates to resonant radiation. Room temper- 
ature spectral narrowing of emission!'! and cooperative emission00 in 7r-conjugated poly- 
mer thin films, superradiant lasing from the J-aggregated cyanine dye molecules adsorbed 
onto colloidal silica and silver^E as well as the room temperature polariton emission from 
strongly coupled organic semiconductor micro cavities^ have been reported. All these effects 
unambiguously mean a collectivization of polymers and J-aggregates via the emission field 



and are of great importance from the viewpoint of laser applications (see Refs. |15Hl7| for 
reviews). On the other hand, the conditions of such a strong coupling are those necessary 
for manifestation of a bistable behavior of a collection of homogeneously broadened two- 
level systems-EHnl Much efforts have been undertaken to this problem with its projection to 
J-aggregated assemblies. The conditions for the bistable optical response to observe from 



an 



individual J-aggreg ateEHi and even of a dimerEM have been analized. Bistability 



discussed in Refs. |2^- p7| , p0| is attributed to an individual aggregate and consists of an abrupt 
switching of the aggregate population from a low level to a higher one as the pump inten- 
sity rises. Consequently, the reflectivity (transmittivity) of a macroscopic ensemble of such 
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species may be switched off (on) in a step-wise way, allowing thus creation of an all-optical 
bistable element. 

The effect we are speaking about originates from the dynamical resonance frequency shift 
dependent on the population of the system. However, as it has been shown in Refs. p5|.[27| 
and [H], an aggregate of a length smaller than the emission wavelength does not display any 
bistable behavior. This finding makes the above mechanism of the aggregate bistability to 
be hardly experimentally verified. Indeed, despite an aggregate may incorporate to itself 
thousands of molecules and thus may have a length larger than the emission wavelength, in 
reality, only a subsystem of segments, consisting of a portion of the aggregate (the so-called 
coherently bound molecules), responds to the action of a resonant external field.0 Because of 
a disorder (of static and thermal nature) of the surroundings, the exciton coherence length, 
being the physical length of an aggregate in the absence of disorder, is reduced to a size 
N* dependent on the degree of disorder. For the real conditions, the number N* is usually 
smaller than the emission wavelength counted in the lattice unit. At room temperature, 
N* is determined basically by dephasing (originated from thermal fluctuations of molecular 
positions) and has an order of magnitude of several lattice units.ilil On decreasing the tem- 
perature, disorder tends to be of static nature. It is determined by static fluctuations of the 
aggregate structure as well as the surroundings and becomes dominating as the temperature 
approaches zero. The static disorder results in Anderson localization of the exciton within 
an aggregate segment with the typical size of the order of several tens of lattice units.fHIBIl 
The number of molecules within the localization length plays now the role of N* and, in fact, 
represents a higher limit for the number of coherently bound molecules. Further increasing 
iV* would be possible by means of reducing the degree of static disorder that is generally 
out of control. 

In this connection, in Refs. ^ the question was risen whether an ensemble of localized 
exciton states may behave in a bistable fashion. It was used the fact that the exciton states, 
located within the same localization domain and being active in the optical response, form a 
local energy structure with a few levels (two or three) similar to the one existing on a regular 



chain with length N ~ iV*.Eill In other words, an aggregated sample can be modelled by an 
ensemble of short (compared to the emission wavelength) linear Frenkel chains stochastically 
distributed over their lengths N. It was shown that an ultrathin film with a thickness much 
smaller than the emission wavelength (which makes possible to use the mean-field approach) 
manifested bistability of the resonance optical response. 



In this paper, we follow the conjecture proposed in Refs. [36] and generalize the results 
reported there to film thicknesses of the order of or larger than the emission wavelength, when 
the mean-field approach, being an adequate approximation for thinner films, is no longer 
valid. This paper is organized as follows. In Sec. [HJ, we present the model and mathematical 
formalism. Section |T| deals with an analysis of the feasibility of a bistable response from an 
ultrathin film (under the condition of validity of the mean field approximation), taking into 
account the fluctuations of all exciton optical parameters, such as the one-exciton transition 
frequency and dipole moment as well as the one-exciton exciton radiative rate.0 Further 
(Sec. [TV] ), we present results of numerical simulations for thin films of thickness larger than 
the emission wavelength and determine the ranges of driving parameters, for which the film 
reflectivity and transmittivity manifests bistability. In Sec. |V|, we make estimates for J- 
aggregates of PIC and show that critical parameters for the occurence of such a behavior 
seems to be achievable for thin films at low temperatures. Finally, Sec. [VI] concludes the 
paper. Some preliminary results of the present study have been reported in our recent 
paper .0 



II. DESCRIPTION OF THE MODEL 

An elementary object of an ensemble we will be dealing with represents a short (with 
a length smaller than the emission wavelength) ordered linear chain consisting of N two- 
level molecules (with TV much larger than unity). Due to the strong intermolecular dipolar 
coupling, the optically active states are the Frenkel exciton states rather than the states of 
individual molecules. In the nearest-neighbour approximation, which we adopt hereinafter, 
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one-dimensional Frenkel excitons appear to be non-interacting fermionsEil so that any state 
with a fixed number of excitons can be constructed as a Slater determinant of one-exciton 

states 

/ 2 \ 1/2 N nkn 
|Jfe) = (_^_) y s in^-|n), k = l,2,....N, (1) 

where |n) is the ket- vector of the excited state of the nth molecule. The energies of the 
one-dimensional exciton gas may take the values W = J2k=i n kEk, where n& = 0, 1 is the 
occupation number of the kth one-exciton state and Ek is the corresponding energy given 
by 

E k = huj 21 - 2U cos -^j- , (2) 

where u>2i is the transition frequency of an isolated molecule and U (chosen hereinafter to 
be positive) is the magnitude of the nearest-neighbour hopping integral. 

The optical transition from the ground state of the chain to the lowest state of the one 
exciton band k = 1 has the dominating oscillator strength (81% of the total one; see, for 
instance, Ref. |6|). Furthermore, the transition betweeen the bottoms of one-exciton and two- 
exciton bands is blue-shifted as compared to the ground-state-to-one-exciton band bottom 
transition by an energy E2 — E\ = 3tt 2 U/N 2 . The blue shift originates of the fermionic 
nature of one-dimensional Frenkel excitons and was experimentally verified for the first 



time in Ref. ^2|. These two facts were put forward in Ref. |36] as a motivation to consider 
the transition from the ground state to the bottom of the one-exciton band as an isolated 
two-level transition. This assumption can indeed be approved provided that the actual Rabi 
frequency of the external field is smaller than the blue shift of one-to-two exciton transitions. 
Further analysis of the transitions to higher exciton manifolds showed that they did not lead 
to qualitative changes of the behavior found within the framework of the two-level model.0 
Based on these findings and taking into account that the spatial inhomogeneity of the field 
and matter variables complicates the treatment drastically, we will use in what follows the 
two-level approach to describe the matter response. 
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Thus, we model the ensemble of linear Frenkel chains as that comprised of inhomoge- 
neously broadened (due to the fluctuations of chain lengths) two-level systems, but with 
the optical characteristics of the transition having all attributes of a one-exciton transition. 
Note that neglecting the one-to-two-exciton optical transitions means that no more than 
one exciton per chain is created by the external field. Respectively, such a channel of effi- 
cient exciton quenching as the intra-chain exciton-exciton annihilation^ plays no role within 
the framework of the two-level approximation. In principle, the excitonic annihilation can 
involve two excitons created on different chains. As follows, however, from the results pre- 



sented in Ref. [44|, the rate of such a process is negligible under the restrictions of our model: 
low temperature and N 3> 1. 

We also assume that the transition dipole moments of all chains are parallel to each 
other as well as to the film plane.@ With regard to the input field £j, the on-resonance and 
normal incidence conditions are chosen. The input field polarization can be set without loss 
of generality to be directed along the transition dipole moment. Then, all quantities can be 
considered as scalars. 

Under the above limitations, the time evolution of the film is described in terms of the 
2x2 density matrix p a p (a, (3 = 1,2) standing for determination of the state of a chain of 
length N. The density matrix equation together with the Maxwell equation for the total 
field £ , including a secondary field produced by the film, form the closed system of equations 
in the problem we are dealing with. It reads 

P2i = -{iu} + T)p2i-i—Z, (3a) 
dEs 

Z = 2t—[ Pl2 - p 21 ] - Y X {Z + 1) , (3b) 
n 



£(x,t) = £Ax,t) / dx'-^-V (x',t 

c Jo at \ 



(3c) 



where the dots denote time derivatives; Z = p 22 — pn; u> = u 2 i — 2{U/K) cos[7r/(iV + 1)] rs 
uj 2 i — 2U/h + Uir 2 /hN 2 is the transition frequency for an individual chain of size iV; d is 



the transition dipole moment of a chain of size N scaled as d = d^yN , where do is the 
transition dipole moment for an isolated molecule; Ti is the spontaneous emission constant 
of the optically active one-exciton state: r\ = jqN with 70 being the analogous constant for 
an isolated molecule (for the sake of simplicity, we have replaced the numerical factor 8/71" 2 
in the expression for d and r\ by unity); r = Ti/2 + T 2 is the dephasing constant including 
the contribution (T 2 ) not connected with the radiative damping. 

The last formula of the set (§) is nothing else but the integral form of the Maxwell 
equation for a film, in which c and L stand for the speed of light and for the film thickness, 
respectively, and V is the electric polarization: 

V = n (d{p 21 + p 12 )) , (4) 

where no is the volume density of chains in the film and the angle brackets denote averaging 
over the chain length distribution with a probability distribution function p(N): (...) = 
EnP(N).... 

As was pointed out in Refs. pi|P^ , |47| -pT|, using the integral wave equation to study 



the non-stationary nonlinear response from a dense absorber has an obvious advantage: 
the boundary conditions are not required within this framework. Indeed, the reflected 
£ r = £(0,t) — £i(0,t) and transmitted £ tr = £(L,t) fields can be calculated obviously on 
the basis of Eq.(|3c|) if the spatial distribution of the electric polarization V is known. The 
latter, in turn, is calculated from constituent equations (|3a|) - ( |3T5| ) and (|J). The background 
refraction index is not included in Eq.fl3c|) because, in fact, this simply results in constant 
renormalization.^1 

Rigorously speaking, for a dense system as that we are considering, the deviation of the 
acting field from the average (Maxwell) field may be significant. A simp least way to account 
for this difference is to add the (Lorentz-Lorenz) local-field correction in the form (47r/3)V 
to the Maxwell field £ in the matter equations (|3a])-(|3b"|) (see, for instance, Ref. |21]). Along 
this paper, however, we will neglect the local field correction justifying this approximation 
by our previous study: for that arrangement of the incident fequency we are going to study 
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(see Sec. |I V B|) , this correction leads only to small quantitative effect s.E! 

Assume that the incident field has the form Si = Ei(t) cos(u){t — kix), where u>i and 
ki = Ui/c are the frequency and wavenumber, respectively, and E{(t) is the amplitude, slowly 
varying in the scale of the optical period lisjiOi. Thus, by passing to the rotating frame by 
means of the representation p 2 i — — (i/2)i?exp(— iu { t) and £ = (l/2)£7exp(— lurf) + c.c, 
with complex amplitudes R and E slowly varying in time, and neglecting the counter-rotating 
terms, we obtain the set of truncated equations 

R = -(iA + T)R + ^Z , (5a) 
Z=- — (ER* + E*R)-V 1 (Z + 1) , (5b) 

E(x, t) = Ei(t- -^j e ikiX + 27m ki £ dx r e ikilx ~ x '\ (dR (x', t - ^ ) , (5c) 

Here, the notation A = uj — cjj is introduced. 

The reflected and transmitted waves we are interested in are simply given by the formulae: 



EJt) = £7(0, t) - E (t) = 2-KkiTiQ [ L dxe 

Jo 



(6a) 



E t {t) = E(L,t) 



Eiit) + 2-KkiUo £ dxe~ ikiX (dR (x,t- 



L — x 



(6b) 



As in our previous papers,^h0 we will neglect the retardation effects replacing t — \x — x'\/c 
by t since the film thickness of our interest is of the order of a few vacuum wavelengths, 
Aj = 2n/ki. Thus, the actual passage time of the light through the film has the order of 
magnitude of the optical period while all the characteristic times of the problem discussed 
(T, h/dEi, A" 1 ) are much longer. One can, therefore, consider the field as propagating 
instantaneously inside the film. 

To perform numerical calculations, let us rewrite Eqs. (|5|) in a dimensionless form using 
the dimensionless field amplitudes e = dE/hT and = dEi/TiT as well as the dimensionless 
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spatial and temporal coordinates £ = kiX and r = ft, where d = (d) and T = (T) are the 
mean transition dipole moment and the mean relaxation constant, respectively. One thus 

gets 



R = -(id + j)R + }teZ 



(7a) 



Z = -~ti{eR* + e*R) - lx {Z + 1) , 



e(£,r) = ei (r)e l « d^e^^Ri^r 



e r (r) = * J"* 1 dZe*(nR(£,T) 



et(r) 



k, L 



ei (r) + * / dte-*UR{t,r) 



(7b) 
(7c) 
(7d) 
(7e) 



where 5 = A/r , 7 = T/r , // = d/d , 71 = iyr, and * = 2^d 2 n Q /hT . 

In our simulations, we choose as distribution function p(iV) a Gaussian centered around 
iV with standard deviation a 



p{N) 



2na 



exp 



(N - N)' 
2a~ 2 



(8) 



assuming a < N. Under this assumption, it is easy to show that the distribution function 
of the detuning S also presents a Gaussian centered at 5q = {uo 2 i — 2U /% + Utt 2 /?iN 2 — u)/T 
with standard deviation a = 2n 2 Ua/hTN 3 . The quantity a can be identified with the 
inhomogeneous width of the exciton absorption line, so that the limits a C 1 and o ^> 1 
(or 2n 2 Ua/hN 3 <C f and 27i 2 Ua/hN 3 3> f in dimensional units) correspond to the cases of 
dominating homogeneous and inhomogeneous broadening, respectively. 

Equations (|?p constitute the basis of our analysis of the optical response from a thin film 
consisting of linear Frenkel chains. We will be interested in the reflection and transmission 



coefficients for amplitude, which are given by TZ 
e(kiL,T)/ei(kiL,T) . 



[e(0,r) -e 4 (0,r)]/e 4 (0,r) 



and T 
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III. ULTRATHIN FILM, L < A' 



We turn first to the case of an ultrathin film (L < A' with A' being the wavelength inside 
the film) when the spatial dependence of R and Z can be neglected and the exponentials in 
Eq. (fTcf) can be replaced by unity inside the film. Then, equations for the fields are simplified 
drastically: 



e(r)=e t (r)+^R(r)) 



(9a) 



e r (r)=V(^(r)) 



(9b) 



e t (r) =e,(r) + ^{fiR{r)) 



(9c) 



where ip = ^>kiL. 

First of all, we are interested in the stationary states in which the system can be found. 
This implies to look for steady-state solutions to Eqs. ([7|), i.e., letting R = Z = 0. Under 
the simplification @, valid for an ultrathin film, it is straightforward to arrive to a closed 
equation for the transmission coefficient T: 

2 



T 



7 



5 2 + 7 2 + /i 2 e?( 7 /7i)^ 2 
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5 2 + 7 2 + /i 2 e 2 ( 7 /7i)^ 2 
Whenever T is found, the reflection coefficient TZ can be expressed through T as follows: 



(10) 



K 2 = V 2 



7 



<5 2 + 7 2 + /i 2 e 2 ( 7 / 7l )T 2 



5 



r 2 



(ii) 



(5 2 + 7 2 + /i 2 e 2 ( 7 / 7 i)T 2 , 

In absence of the chain length fluctuations, Eq. ( |TD| ) is of third order with respect to T 2 and 
thus may have three real roots,i! indicating the possibility of bistable behavior of the system 
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transmittivity. It occurs sX ip > ip c = Due to the relation (|TT|), the system reflectivity 
is expected to behave in the same manner. The chain length fluctuations make in general 
impossible to predict the order of the resulting (after carrying out the averaging procedure) 
transcendent equation. 

In Ref. [55], Eq. ( |IT)D was analyzed under the assumption that the main effect of the 
length fluctuations originates in the distribution of detuning 5, while fluctuations of the rest 
of parameters (/i, 7, 71) are of no importance. Replacing then the distribution over length, 
p{N), by a distribution over detuning and taking the latter in the form of a Lorentzian 



of width g, authors of Ref. |36| could evaluate integrals in Eq. ( |10D analytically and found 
that the resulting equation gaves rise to a bistable behavior, even under the condition of 
dominating inhomogeneous broadening [g > 7). 

Here, we do not restrict ourselves to the above assumption and carry out calculations, 
accounting for fluctuations of all stochastic variables in Eq. ([101) . We present the results 
for a particular case of the incident frequency tuned up to the center of the absorption 
band (#o = 0). Figure [TJ shows a "phase" diagram of the system behavior in the space of 
parameters (a, ip) calculated on the basis of Eq. flT0|). The parameters of Gaussian p(N) 
were chosen as N = 30 and a = 9. By the terms "one solution" and "three solutions" in 
Fig. [l], we marked regions where Eq. fllQl) had single-valued and three-valued real solutions, 
respectively. The solid curve, which separates these two regions, is nothing but the critical 
value of ip for bistability to occur versus the inhomogeneous width o = 2ir 2 Ua/hTN 3 . It 
starts from tp c = 8 at a = 0, in full correspondence with the earlier findings, and then 

gradually goes up on increasing a. 

For comparison, we also depicted in Fig. [1] the curve (dashed) obtained under the ap- 
proximations used in Ref. |36j, i.e., assuming the detuning 5 as being the unique stochastic 
variable and setting /1 = 7 = 1,71 = const. As can be seen, there is almost no difference 
between the solid and dashed curves. This means that indeed, fluctuations of the detun- 
ing basically determine the result of averaging in Eq. flTU|), thus approving this heuristic 



assumption used in Ref. [36 
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Figure presents the typical examples of the input field dependence of the reflection 
and transmission coefficients for an ultrathin film {h^L = 0.1) which were obtained by the 
numerical solution of Eqs. ([/]) at adiabatic scanning of the input field amplitude up and 
down. The results were obtained chosing the following set of parameters: ip = 20, N = 30, 
a = 9, 71 = 0.25N/N, 7 = 0.875 + 0.125N/N. From this figure, one can conclude that, 
when the inhomogeneous width a is small, the system shows a stable hysteresis loop (optical 
hysteresis) both in reflectivity and transmittivity, despite the fact that fluctuations of all 



parameters are taken into account (note that the hysteresis of transmittivity in Ref. |36 



was 



calculated only for the particular case of absence of the chain length fluctuations). However, 
as a grows, the hysteresis cycle disappears and only a one-valued monotonic response exists. 

In order to gain insight into the time required for the output signal to approach its 
stationary value, we have numerically solved Eqs. (0) at a fixed amplitude of the input field, 
ej, and its further step- wise switching to another value. Figure |3] shows an example of such 
calculations for the same set of parameters as in Fig. |2|, setting a = 0.25. We start with 
a non-saturating incident field amplitude = 3, for which the reflection (transmission) 
coefficient is high (low), and wait for a stationary value of the latter. At some instant of 
time (in particular, at r = Tt = 50), we switch step- wisely the incident field to a saturating 
value Ci = 6, for which the reflection (transmission) coefficient is low (high), and wait again 
for a stationary value of the latter. The results of these calculation are depicted in Fig. |3| 
with solid lines. With dotted lines, we present the results of analogous calculations obtained 
only by switching back the incident field amplitude from a higher (e$ = 6) to a lower (e^ = 3) 
value. 

From these data, it can be claimed first that the transient time depends on whether the 
incident field amplitude is switched up or down. Second, the transient time is of the order of 
a few units of the population relaxation time T^ 1 . Indeed, in units of f _1 this time interval 
is of the order of 10. Taking into account that T\/T = 0.25 in this particular calculation, 
one arrives at the above stated conclusion. 
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IV. THICK FILM, L > X' 



A. Motivation 

In our previous studies of nonlinear resonant reflection from an extended {L > A') dense 
system of homogeneously broadened two-level molecules,&S we found that the reflectivity 
of such a system could show different regimes (bistability and self-oscillationsE§E§ as well as 
chaos^) basically governed by the parameter instead of ij) = ^kL for an ultrathin film 
(see also Refs. |52"| , [53| ). The physical meaning of \l/ is nothing but the halfwidth (in units of the 
homogeneous width F) of a gap in the spectrum of field + molecules collective excitations - 
polaritons (polariton splitting; see for more details the monograph by DavydovCT). This gap 



appears in the vicinity of the molecule-field resonance and requires the condition ^ ^> 1 to 
be fulfilled. The linear dielectric constant for frequencies ranging within the gap is negative 
(the refractive index is imaginary), implying total reflection of the light of such frequencies.il 
The physical origin of the effects reported in Refs. |48|- |50| , p3| is attributed to saturation 
of the refraction index by the field acting within the polariton band gap. Let us recall 
briefly the motivation rised in those papers. For an input field varying slowly in the scale 
of the relaxation times T -1 , Tj" 1 (the case of our interest in the present study), the medium 
adiabatically follows the field (R can be set to zero). Under such conditions, Eqs. ( |7a]) 
and (|?c| ) are equivalent to only one equationJiH 



2 +e(|e| 2 )e = , (12a) 



The quantity £(|e| 2 ) is the field-dependent dielectric function in which fluctuations of all 
parameters (/x, 7,71,5) are explicitly taken into account. The limit of absence of the chain 
length fluctuations (/x = 7 = 1, 71 = const and 5 = const) formally corresponds to the case 
considered in Refs. 
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£(N ^ 1+2 w i+ ;^ e|2/7i . d3) 

We assume now that \l/ S> 1 and neglect for a moment the field (linear case). Then, as can be 
seen from Eq. (|13D , the linear dielectric function e has a dominating negative real part within 
a rather wide interval of changing the detuning, in fact, from a few negative units to — 2\&. 
A weak field acting within this band will be totally reflected from the system boundary 
At higher amplitudes of the field, 7|e| 2 /7i ^> the dielectric function is saturated by the 
acting field, i.e., goes to unity, producing conditions for penetration of the field into the 
medium. The higher the field amplitudes, the deeper the field penetrates into the medium. 
Now, reflection occurs from a very narrow (due to a high magnitude of \&) interface between 
the saturated and non-saturated regions.^! This interface plays the role of an "artificial" 
mirror giving rise to a feedback necessary for the effects outlined above to build up. 

It turned out that namely in the limit of larger i.e., in the presence of the polariton 
band gap, the system shows the above mentioned peculiarities of nonlinear optical response 
that are absent in the opposite case \P < 1. When inhomogeneous broadening dominates 
(the main range of our interest), one should compare the inhomogeneous width 2cr with the 
width of the gap 2\I/. A simple analysis of Eq. (|12b ) reveals the clear result that a well- 



pronounced gap (where the dielectric function has a dominating negative real part) develops 
at \I/ 3> o. An inhomogeneous broadening of the order of, or higher than 2\1/ destroys 
the gap and, as a consequence, all the nonlinear effects accompanying the presence of this 
gap. Below, we show details of the influence of inhomogeneous broadening on the nonlinear 
optical response of the film. 

B. Numerical simulations and discussions 

All nonlinear features accompanying the presence of the polariton gap were found to be 
stronger when the input field frequency lies in the middle of the gap, i.e., at 5 = — vl/.&ilH 
Therefore, when studying the nonlinear reflection from and transmission through a film 
consisting of inhomogeneously broadened two-level systems, we will also set 5 = — \I/ and 
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consider \& > \& c = 4.66 (\l/ c is the threshold for bistability in reflection^) to get the optimal 
conditions for the effects we are looking for. It is worth noting that, as was found in Ref. |49|, 
the local-field correction at 5 = — \I/ has a relatively small quantitative effect. In particular, 
the threshold value \I> C changes from 4.66 to 5.45 if the local-field correction is taken into 
account. 

Figure f| presents examples of the input field dependence of the reflection and trans- 
mission coefficients for a thin film of thickness L — Aj depending on the inhomogeneous 
width a = (2ir 2 a/N 3 ) ■ (U/KT), where N and a were fixed (N = 30, a = 9) and U/ht 
was varied. The other parameters were chosen as follows: \& = —5 = 6, 71 = 0.25iV/iV, 
7 = 0.875 + 0.125N/N. The data were obtained at adiabatic scanning of the input field 
amplitude up and down. From this figure, one can conclude, first, that the system shows 
a stable hysteresis loop both in reflectivity and in transmittivity or, in other words, be- 
haves in a bistable fashion until a < i.e., until the inhomogeneous width approaches 
the polariton splitting, in full correspondence with what we mentioned above. The second 
observation, evident from Fig. f|, is that the switching amplitudes of the input field do not 
depend strongly on the inhomogeneous width (at a fixed magnitude of 

In Fig. |5], we depicted spatial profiles of the amplitude module of the field inside the 
film versus the input field amplitude e«. The calculations were performed for the same set of 
parameters and conditions as those presented in Fig. |], only setting a = 2. Plots (a) and (b) 
correspond to adiabatic scanning of up and down, respectively. The darkness of a local 
differential domain is proportional to the field amplitude module. Observing these plots, we 
can claim the following. Despite the fact that the thickness of the system does not exceed 
one wavelength in vacuum, A«, the field inside the film is not uniform at any amplitude of 
the input field. For not-saturating magnitudes of e$, this is simply explained by the fact 
that the (linear) refraction index in this case (at 5 = n w —1, that, in turn, gives rise 
to decreasing the field amplitude on a scale short compared to Aj. At higher magnitudes of 
the input field, ranging within the bistable interval (see Fig. |]), the spatial inhomogeneity of 
the field demonstrates a global character, originating in the interplay between the forward 
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and backward (in the present case, reflected from the back boundary of the film) waves. 

Comparison of the field profiles depicted in the plots (a) and (b) at a fixed magnitude of 
the input field reveals that they differ one from the other. In other words, one may make a 
statement about the existence of spatial hysteresis of the field inside the film. 

Figure ^ sheds light on the transient time of approaching the reflectivity (a) and trans- 
mittivity (b) their stationary values when switching the incident field amplitude up (solid 
lines) and down (dotted lines), similar to that presented in Fig. |||. The calculations were 
done for the following parameters: ^ = —5 = 6, N = 30, a = 9, a = 2, 71 = 0.25N/N, 
7 = 0.875 + 0.125iV/iV As in the case of an ultrathin film, the transient time has the order 
of several units of the population relaxation time F^ 1 . 

In Fig. [F] we present the data of similar calculations as depicted in Fig. |] performed 
for thicker films: L = 3Aj and L = 5Aj. As follows from these results, the hysteresis in 
reflection remains to be very pronounced. However, it tends to disappear in transmittivity. 
It is simply due to the fact that in order to make transparent a film of higher thickness, one 
needs to apply an input field of very high amplitude, for which the system already does not 
manifest at all bistability in transmission. By contrast, with regards to reflection, the film 
behaves as a " semi- infinite" medium, for which bistability in reflection may occur provided 
the parameter \l/ exceeds its critical value \l/ c = 4.66.@. In Fig. || we depicted the spatial 
profile of the field inside the film with thickness L = 3Aj for the conditions for bistability to 
occur. It is clear from these pictures that the field is absolutely attenuated within a depth 
a bit longer than Aj and, hence, the transmittivity in this case is negligible. 

As follows from our previos study&"0 of the nonlinear reflection from a collection of dense 
homogeneosly broadened two-level atoms, the system may show instabilities of different types 
(self-oscillation and chaos) at higher values of the polariton splitting ty. Figure ^ presents 
an example of such a behavior of the film reflectivity (self-oscillation) for \I/ = 10 and the 
inhomogeneous broadening a = 2. The attempt to get a hysteresis loop of the reflection 
coefficient reveals that the lower branch is unstable (see Fig. [5]a). The calculation done at 
a fixed amplitude of the input field (e* = 8.5, above the switching threshold) shows that 
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indeed, the system has a regime of self-oscillations with a frequency of the order of Y/2. 

V. ESTIMATES OF PARAMETERS 

First, let us discuss the applicability of the two- level model used in this Paper. One 
should compare the typical energy spacing between the first and second exciton levels (E2 — 
Ei)/hT ~ 3n 2 (U/hT)/N 2 with the switching magnitude of the input field ej. For the 
parameters at hand, the former quantity ranges within the interval [3.3, 33] while the latter 
is approximately equal to 3. Therefore, we can conclude that our approach is correct and 
there is no necessity to include the one-to-two exciton transitions. 

It is worthwhile to analyze the parameters of real systems in order to get insight into 
feasibility of the bistable mechanism we are dealing with. In this sense, the J-aggregates of 
PIC, as one of the most studied species of the type we need, seem to be suitable objects. 
The width of the red J-band of PIC-Br (centered at A = 576.1 nm), at low temperature, 
has an inhomogeneous nature and an order of magnitude of 30cm _1 (or lps' 1 in frequency 
units)i~@. The blue shift of the transition from one-to-two exciton bands with respect to that 
from the ground state to the one-exciton band for the red J-band has the same order. The 
polariton splitting is \&r = 2ird 2 n /h = (3/ 16n 2 )j Nn X 3 , so that taking 70 = (l/3.7)ns _1 
and Ntiq = 10 18 cm -3 (an achievable concentration of molecules before aggregation), we 
obtain ~ lps~ l , i.e., a value that exceeds twice the halfwidth of the J-band 0.5ps~~ l (oT 
in our notation). Hence, for the parameters used, we are under the conditions needed for 
bistability to occur. 

VI. CONCLUSIONS AND REMARKS 

In this Paper, we have numerically studied the optical bistable response of a thin film 
with thickness of the order of and larger than the emission wavelength, comprised of oriented 
linear Frenkel chains. We have taken into account a distribution of chains over length, 



18 



resulting in an inhomogeneous broadening of the exciton optical transition as well as in the 
distribution of the transition dipole moment. 

From our study, the following conclusions can be drawn: 

(i) Within a certain range of driving parameters (resonance detuning, polariton splitting, 
and inhomogeneous width) , the films of thickness of the order of the vacuum emission wave- 
length really show a bistable behavior with respect to both the reflection and transmission 
of the resonant light; 

(ii) The bistability effect exists until the inhomogeneous width of the transition ap- 
proaches the polariton splitting, i.e., for a fairly wide range of widths; 

(iii) With increasing the film thichness, the bistability in the transmitted signal tends to 
disappear while it remains in the reflected signal; 

(iv) At higher magnitude of the linear refraction index, the film reflectivity shows insta- 
bilities wich are of self-oscillating type; 

(v) The parameters needed for the bistable behavior seem to be achievable for a thin film 
of J-aggregates of polymethine dyes and for some classes of conjugated polymers deposited 
onto a dielectric substrate such as poly(p-phenylene-vinylene) derivatives. 

The inhomogeneous broadening acts as a destructive factor. Reducing any source of 
inhomogeneity is thus of great importance. In this sense, thin films comprised of thiophene 
oligomers seem to be very attractive objects, too, due to the possibility of precise controlling 
the sizes of this type of molecules.i!~§! 

To conclude, it is to be noting that the authors of Ref. |59] reported a room-temperature 
formation of polariton states in ordered cyanine dye films. The latter thus can also be 
considered as a promising object from the viewpoint of our findings. 
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FIGURES 

FIG. 1. (a, tp)-m&p obtained by solving Eq. (|lO|). The solid line separates two regions within 



which there are one and three real solutions of Eq. (10) and thus represents the critical value of 
ip for bistability to occur versus the inhomogeneous width a = 2it 2 Ua/hT iV 3 . The dashed line is 
calculated under the assumption that 5 is the only stochastic parameter, as was done in Ref. [36]. 

FIG. 2. Amplitude reflection (a) and transmission (b) coefficients of an ultrathin film 
{kiL = 0.1) calculated using Eqs. ([?]) with the input field amplitude ej = dEi/hT scanned up 
and down for different values of a = 2ir 2 U a/hTN 3 where U/hT was varied. Other parameters were 
chosen as follows: tp = 20, N = 30, a = 9, 71 = 0.25iV/iV, 7 = 0.875 + 0.125N/N. Time is in units 

of r _1 . 

FIG. 3. Kinetics of approaching the reflection (a) and transmission (b) coefficients their 
stationary values after a sudden switching (at an instant r = Tt = 50) of the incident field 
amplitude. The calculations were done for the parameters of Fig. Q setting a = 0.25. The solid 
lines represent the results obtained when the input field amplitude was switched from a low value 
ei = dEi/hT = 3 (a high reflection) to a higher value ej = dEi/hT = 6 (a low reflection). The 
dotted lines represent the results for the back switching of the input field amplitude, from = 6 
to e« = 3. Time is in units of F _1 . 

FIG. 4. The amplitude reflection (1Z) and transmission (T) coefficients of a film with thickness 
L = Xi (kiL = 2tt) calculated at adiabatic scanning of the input field amplitude = dEi/hT up 
and down for different values of the inhomogeneous width a = 2ir 2 Ua/hTN 3 , where U/hT was 
varied. Other parameters were chosen as follows: * = 2ircPn /hT = -S = 6, 71 = 0.25iV/iV, 
7 = 0.875 + 0.125iV/iV. Time is in units of f _1 . 

FIG. 5. Examples of the spatial profiles of the field amplitude module inside the film corre- 
sponding to the case with a = 2 in Fig. [|. The graphs (a) and (b) show the inside field profiles when 
the input field amplitude ej = dEi/hT is scanned up and down, respectively. The darkness of a 
local differential domain is proportional to the module of the inside field amplitude, |e| = d\E\/hT. 
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FIG. 6. Kinetics of approaching the reflection (a) and transmission (b) coefficients their 
stationary values, calculated after a sudden switching (at an instant r = Tt = 50) of the input 
field amplitude from a value ej = dEi/fiT = 2 below the switching point ~ 2.9 to a one 
ei = dEi/fiT = 3.5 above the latter (solid line). By the dotted line, the back switching is shown. 
The parameters are the same as in Fig. [|. Time is in units of f 1-1 . 

FIG. 7. The same as in Fig. || (a) for the case of a = 2 and several values of the slab thickness. 

FIG. 8. The same as in Fig. |5| for a slab of thickness L = 3Aj. 

FIG. 9. a - Optical hysteresis loop of the reflection coefficient for a film of thickness L = 2Aj 
revealing an instability. The calculation was done by means of Eqs. (|?|) at adiabatic scanning of the 
input field amplitude e, = dEi/fiT up and down, setting the following parameters: ^ = 10 = —5, 
N = 30, a = 9, a = 2, 71 = 0.75N/N, 7 = 0.625 + 0.375(iV/iV). b - Kinetics of the reflection 
coefficient 1Z calculated at a fixed amplitude of the input field (e^ = 8.5) showing self-oscillations. 
The set of parameters is the same as in plot a. Time is in units of r . 
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